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1 Introduction 

This worksheet demonstrates the use of Maple in Linear Algebra. 

We give a new procedure (PowerMatrix) in Maple for finding the fc th power of n-by-n square 
matrix A, in a symbolic form, for any positive integer k 7 k > n. The algorithm is based on an 
application of Cayley-Hamilton theorem. We used the fact that the entries of the matrix A fc satisfy 
the same recurrence relation which is determined by the characteristic polynomial of the matrix A 
(see [1]). The order of these recurrences is n — d, where d is the lowest degree of the characteristic 
polynomial of the matrix A. 

For non-singular matrices the procedure can be extended for k not only a positive integer. 

2 Initialization 

> restart: 

with (Linear Algebra) : 



2.1 Procedure Definition 



2.1.1 PowerMatrix 

Input data are a square matrix A and a parameter k. Elements of the matrix A can be numbers 
and/or parameters. The parameter k can take numeric value or be a symbol. The output data is 
the fc th power of the matrix. The procedure PowerMatrix is as powerful as the procedure rsolve. 
> PowerMatrix := proc(A: :Matrix,k) 

local i, j ,m,r,q,n,d,f ,P,F,C; 

P := x->CharacteristicPolynomial (A,x) ; 

n := degree (P (x) , x) ; 

d := ldegree(P(x) ,x) ; 
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F := (i , j ) ->r solve ( sum (coef f (P(x) ,x,m)*f (m+q) ,m=0. .n)=0,seq(f (r)=(A~r) [i , j] , 

r=d+l. .n) ,f) ; 

C := q->Matrix(n,n,F) ; 

if (type (k, integer) ) then return(simplif y (A~k) ) elif (Determinant (A) =0 and 
not type (k, numeric)) then printfO'The 7,a th power of the matrix for %a>=°/.d: ", 
k,k,n) elif (Determinant (A) =0 and type(k, numeric)) then return (simplif y (A~k) ) fi; 
return ( simplif y( subs (q=k,C(q)) ) ) ; 
end: 



3 Examples 



3.1 Example 1. 

> A := Matrix ([[4,-2,2] , [-5,7,-5] , [-6,6,-4]]); 



A:= 



4 
-5 
-6 



2 

-5 
-4 



> PowerMatrix(A,k) ; 



-2 fc + 2-3 fe 
-5-3 k + 5-2 k 
6 • 2 k - 6 • 3 fc 



> Determinant (A) ; 



> B:=A~(-1); 



2 (i+fe) -2-3*= 
5 • 3 fe - 4 • 2 k 
-6 • 2 fe + 6 • 3 fe 
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_ 2 (i+fe) +2-3*= 
-5-3 fe + 5-2 fe 
-6-3 fe + 7-2 fe 



B:= 



> PowerMatr ix (B , k) ; 



-2(~ fe ) + 2 • 3 ( ~ fc) 2( 1_fe ) - 2 • 3 ( ~ fe) 
-5 • 3(- fc ) + 5 • 2<- fc ) 5 • 3^ fc ) - 4 • 2^ 
-6 • 3^ fe ) + 6 • 2<- fc ) -6 • 2(~ fc ) + 6 • 3^ fc ' 



-2( 1 - fe ) + 2-3(- fc ) 
-5 • 3^ fe ) + 5 • 2<~ fc ) 
-6 • 3(- fe ) + 7 • 2<~ fc ) 



3.2 Example 2. 

> A:= Matrix([[l-p,p] , [p,l-p]]) ; 



1 ~P P 
P 1 ~P 
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> PowerMatrix(A,k) ; 



(l-2p) fc 1 (l-2p) fc 1 ' 
2 2 2 2 



(l-2p) k 1 (l-2p) k 1 



2 2 2 2 J 



The example is from [4], page 272, exercise 19. 



3.3 Example 3. 

> A:= Matrix([[a,b,c] , [d,e,f] , [g,h,i]]) ; 



A:= 



a b 

d e 



f 



c 



9 h 



> PowerMatrix(A,k) [1,1] ; 

S_R = RootOf ( (gbf + hdc + iea - gee - hfa - idb)_Z 3 + (gc + hf + db - ie - ia - ea)_Z 2 + (i + e + a)_Z - 1 ) 



-3_R 2 3 ce - 3_R 2 hfa - 3_R 2 idb + 2_Rgc + 2_Rhf + 2_Rdb - 2_Rie - 2_Ria - 2_Rea + i + e + a)_R) 

# Warning! 

In this example MatrixPower and MatrixFuction procedures cannot be done in real-time. 

# MatrixPower (A, k) [1,1] ; 

# MatrixFunction(A,v"k,v) [1 , 1] ; 

3.4 Example 4. 

> A:=Matrix([[0,0,l,0,l] , [1 ,0 ,0, , 1] , [0 ,0, , 1 , 1] , [0 , 1 ,0, , 1] , [1 , 1 , 1 , 1 ,0] ] ) ; 




A:= 




1 



1 







1 
1 



1 





1 



1 

1 

1 1 

1 

1 



> PowerMatrix(A,k) [1,5] ; 




Replace ' : ' with ' ; ' and see result ! 

> MatrixPower(A,k) [1,5] : 

> assume (m: : integer) : simplify (MatrixPower (A, k) [1,5]): 
The example is from [3] , page 101 . 
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3.5 Example 5. and Example 6. 

Pay attention what happens for singular matrices. 
3.5.1 Example 5. 

> A:=Matrix([[0,2,l,3] , [0,0,-2,4] , [0,0,0,5] , [0,0,0,0]] ) ; 



A:= 






2 


1 


3 








-2 


4 











5 















> PowerMatrix(A,2) ; 









-4 


13 











-10 



























> PowerMatrix (A , 3) ; 












-20 







































> PowerMatrix(A,k) ; 

The k th power of the matrix A for k > 4: 

" " 







> MatrixPower (A,k) ; 

Error, (in LinearAlgebra: -LA_Main: -MatrixPower) 
power k is not defined for this Matrix 

> MatrixFunction(A, v~k,v) ; 

Error, (in LinearAlgebra :-LA_Main:-MatrixFunction) 
Matrix function v fe is not defined for this Matrix 

The example is from [2], page 151, exercise 23. 



3.5.2 Example 6. 

> A:=Matrix([[l,l,l,0] , [1 , 1 , 1 , -1] , [0, , -1 , 1] , [0 ,0 , 1 , -1] ] ) ; 



A:= 



1 1 

1 1 







1 

1 -1 

-1 1 

1 -1 
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> PowerMatrix(A.k) ; 

> The fc th power of the matrix for fc>4: 



2 (-i+fc) 2 (-l+fe) 

2 (-i+fc) 2 (-i+fc) 






(_l)(l+fe).2 fc 5-2 fc 

1 

16 16 

5-2 fe 5 ■ (-l)( 1+fe ) ■ 2 fe 5-(-l) fe -2 fc 2 fe 
16 + 16 16 16 

(— l) fe • 2 ( ~ 1+fc ) 



2 fe (_i)fc.2 fc 
16 + 16 



16 

(_l)(i+fe) . 2(- 1+fc ) 
(_!)fc . 2 (-i+fc) 



> MatrixPower (A,k) ; 

Error, (in LinearAlgebra:-LA_Main: -MatrixPower) 
power k is not defined for this Matrix 

> MatrixFunction(A, v~k,v) ; 

Error, (in LinearAlgebra:-LA_Main:-MatrixFunction) 
Matrix function v k is not defined for this Matrix 
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5 Conclusions 

This procedure has an educational character. It is an interesting demonstration for finding the k th 
power of a matrix in a symbolic form. Sometimes, it gives solutions in the better form than the 
existing procedure MatrixPower (see example 4.). See also example 5. and example 6., where we 
consider singular matrices. In these cases the procedure MatrixPower does not give a solution. The 
procedure PowerMatrix calculates the fc th power of any singular matrices. In some examples it is 
possible to get a solution in the better form with using the procedure allvalues (see example 3.). 
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